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I. INTRODUCTION 

Since the discovery of Hawking radiation lU in the early 
seventies much research into quantum field theory on gen- 
eral classical backgrounds has taken place. One of the most 
well established results of this field is that of the Unruh effect 
10. HI 01 which states that a constantly accelerated particle de- 
tector in the Minkowski vacuum will register a thermal bath 
of particles at a temperature given by; 



where a is the acceleration of the detector 

Roughly one decade after the Hawking result, Bell and 
Leinaas took these observations and applied them to linearly 
and circularly moving electrons [5] addressing the anomalous 
spin depolarization that was observed in storage rings (see 
also H). 

Mainly due to the drive for quantum computing, the last 
twenty years has also seen the rise of quantum information 
theory, at the heart of which lies the resource of entanglement. 
Entanglement is a phenomenon unique to quantum systems. 
Two systems are said to be entangled if they are described 
by a single state vector that cannot be written as a product of 
state vectors for each part. The entanglement between two 
qubits can be quantified by a function called concurrence |(7(], 
taking values on the interval [0,1], where a concurrence of 
zero (one) represents no (maximal) entanglement. 

It has only recently been shown that the resource of 
entanglement is frame dependent |8], in the sense that the 
distribution of spin and momentum entanglement between 
two electrons can depend on the speed of the observer 
measuring these properties. Furthermore, it was shown by 
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llH [T^ that the amount of entanglement shared by a pair of 
field modes in a constantly accelerating cavity entangled with 
another pair at rest decreases with acceleration. For further 
developments in this field see llTll[T2l[T3ll . 

While the entanglement between electrons under Lorentz 
transformations and the entanglement between Dirac field 
modes in accelerating cavities have both been investigated, 
a systematic study on the spin entanglement between an 
accelerating electron and one at rest remains to be done. At 
first sight there are two effects that require consideration: 

The first effect is a rotation of the accelerating electron's 
spin as viewed from the stationary electron frame. Uniform 
acceleration consists of a continuous sequence of infinitesi- 
mal boosts, dv' = adr, where the acceleration a and change 
in velocity dv' are measured in the Instantaneous Rest Frame 
(IRF), i.e., the frame that is momentarily at rest with the mov- 
ing electron at proper time t. Viewed from the unaccelerated 
frame the change in motion comprises of a Lorentz boost 
A{v) into the IRF followed by an infinitesimal boost A{dv'). 
It is well known that if two successive non-collinear Lorentz 
boosts take X to X' then the single Lorentz transformation 
relating these frames is not a pure boost, but rather is the 
product of a boost and a rotation. The rotation induces an 
operation on the quantum state called a Wigner rotation. As 
discussed in I's', T?] when the particle is not in a momentum 
eigenstate these Wigner rotations transform the entanglement 
between the spin and momentum degrees of freedom. 

The second effect arises from interactions with the thermal 
environment ([T]). As shown in |5], the spin of the accelerat- 
ing electron will flip in response to the Rindler radiation. To 
understand what effect this has on the spin entanglement it is 
instructive to consider the simpler spin flip interaction: 

H,nt--b^i){U + b\^){i\, (2) 

where 6^^' are the annihilation (creation) operators for the 
field modes of the thermal field and the Hamiltonian is only 
acting on the subspace of the first spinor. A first order inter- 
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action will turn the entangled state: 

l*> = 7l(ITT) + IU»|no), (3) 

into: 

l*)^7i(Ut)K + i) + in)K-i)). (4) 

Since one is only observing the spin degrees of freedom the 
field subspace must be traced out. One then finds the spins 
disentangled, i.e.. 



TiXTi 



itXit 



(5) 



In this paper we show that both of these effects can be in- 
corporated into a single framework by using an open quantum 
system formalism and extending the single electron work of 
Bell and Leinaas to an entangled two particle system. In our 
case we take an accelerating electron similar to the one in jst], 
and entangle it with an ancillary rest particle. We will show 
that the first effect (described above) is nothing other than the 
Thomas precession of the electron, which is automatically 
present as a term in the effective Hamiltonian [5]. In section 
Uni we show that the first effect can be neglected for the 
classical path under our consideration. The second effect is 
then the dominant source of disentanglement and within the 
open system framework we show that analytical values for 
the relaxation timescales and the concurrence of the system 
can be obtained. 

It is found that due to the Unruh radiation that the acceler- 
ated particle experiences, the acceleration increases the rate 
at which the entanglement is destroyed. In the limit when 
the Rindler temperature ^ is larger than the spin energy 
separation scale this dependence takes on a particularly 
simple form. For an initially maximally entangled system 
the time to completely disentangle, as measured by the 
accelerated spinor, is found to be: 



accelerations, a(r). It is for this reason we expect that the 
method we use to calculate concurrence in this paper will 
also be applicable to more general accelerations, for instance 
sinusoidal motion. 

Consider a particle confined to move along the z-direction 
and subjected to time-dependent accelerations in its IRF . In 
that frame the particle will gain a small velocity dv' = a{T)dT 
in a small time interval dr due to its acceleration, where a(T) 
is the acceleration as measured in the IRF. A static observer 
will measure a velocity 



dv 



dv' 



1 + vdv'/c^ 



(8) 



where v is the velocity of the IRF at time t(r). Keeping terms 
to order 0{dv') one finds dv = {1 — {v/cY)dv' and that the 
rapidity is given by: 



r(r) = arctanh 



1 



v{t)_ 
c 

dv 



cj 1-(w(t)/c)2 

The particle's trajectory in the static frame is then 

f dz dt 
z(t) = / ——dr, 
^ ' . dtdr ' 



(9) 

a(r)dr, (10) 



(11) 



•j sinh r(T)(iT. (12) 



We can also find the dependence of coordinate time on the 
proper time, ^(t), in flat spacetime using the Minkowski line 
element l = (|f)^-(^)^- This impHes that 



(13) 



^(t) = / cosh r(T)(ir. 



and in the rest frame: 



37rln3 hc^ 
8 ^2^3 ' 



c / Stt In 3 he' 
exp 



2a 



(6) 



(7) 



This paper is organised as follows: after briefly reviewing 
accelerated world lines on flat spacetimes in sectionlUwe out- 
line our configuration of fields and electrons and then deter- 
mine and solve the master equation for this system in section 
Unl In section |IV] we will calculate the relaxation times and 
study the time decay of the concurrence in the system before 
summarising our findings. 



II. ACCELERATED PATHS IN ONE DIMENSION 

In what follows we give a brief derivation of the Rindler 
worldline emphasising its generalisation to time-dependent 



For the special case where the acceleration is constant, a(r) = 
a, we obtain: 



/ \ c a 
t[T) = — smli — T, 
a c 

z{t) — — cosh — T, 
a c 



(14) 
(15) 



which is the usual Rindler result. The path of the Rindler 
worldline is shown in figure[T](a). 



III. QUANTUM EVOLUTION OF SPIN ENTANGLEMENT 
BETWEEN FRAMES 

In our setup, we take two electrons keeping one of them 
(62) stationary and isolated from any fields whilst the other 
(ei) is placed under a constant magnetic field and accelerated 
by a constant electric field, see figure [T] In this situation all 
the dynamics occurs on ei with 62 playing a spectator role. 
62 acts as an entangling partner for ei and also defines a static 
inertial frame. We investigate what effect the acceleration has 
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FIG. 1 : Two different views of the two entangled spins in our setup. 
In a) we have the Minkowski diagram with electron ei following 
the Rindler worldline for a constantly accelerated path while elec- 
tron 62 remains at rest (following the vertical dotted line) and defines 
the inertial frame of this diagram. Figure b) on the other hand is a 
physical view of the two electrons; the top electron, ei, is acceler- 
ated by a constant electric field that points in the negative z-direction 
and placed under a constant magnetic field that points in the positive 
z-direction, while the bottom spinor, 62, is isolated and at rest. 



on any entanglement the two electron spins initially possess. 

Let the electron mass be m, charge be e and the mag- 
nitude of the magnetic moment be /i, and define the fields 
to point along the z-direction i.e., E = (0,0, —i?^) and 
B = {0,0, Bz). The IRF acceleration of the first electron 
ei is given by: 



(16) 



Since the electron charge is negative the electron accelerates 
in the positive z-direction shown in figure[T](b). 

The evolution of the electron spin is determined using clas- 
sical relativistic electrodynamics (see Appendix A and also 



m [Tsl [T^). The effective spin-field interaction Hamiltonian 
for the spin of a relativistic electron 1 17] is: 



H 



h 72 dj3 



X /3 



(17) 



where (5 = v/c, ^ — (1 — fP) ^Z^, cr is the vector of Pauli 
matrices, and B' is the magnetic field in the IRF given by: 



B' = -i\B-l3y.E- 



7 + 1 



(3(3 B 



(18) 



As discussed in Appendix A the second term in equation (fTTb 
is related to the phenomenon of Thomas precession. 

We first analyse this equation semi-classically, treating the 
magnetic field as a classical field and taking the electron path 
to be the classical trajectory given by equations (fT4l) and dTsl l. 
At this level the Thomas term vanishes as the cross product 
between the velocity of the classical particle and its accelera- 
tion is zero. Changing the time parameter to the proper time, 
T, (corresponding to multiplying the Hamiltonian ( fTTI l by 7) 
the evolution of the spin is determined by the usual semi clas- 
sical Hamiltonian 



Hue = -fidzB^ 



(19) 



where we have used the result that B'^ — B^ for the classical 
path and field. This time-independent system has two 
energy eigenstates defined by az = ±1, with an energy gap 
A = 2nBz and therefore spin flipping does not occur at this 
level. 

From the Lorentz force equation (assuming that all non 
electromagnetic forces are negligible), one finds: 



P 



mc{i(3 + -fl3), 
e{E + (3x B) 



and thus equation ( fTTj i can be written as: 

H = -a- I7- VB' + (7 + X 



using: 



E'^-i[E + [3xB- 



7 



-(3{(3 ■ E) 



(20) 
(21) 



(22) 



(23) 



7 + 1 

In terms of the fields in the stationary frame, 

H = -^i<T-{-i-^B + {^+\)-^Ex (3) , (24) 

where /i = 27^^ with g ^ 2, which differs from (|22] | by a 
minus sign in the second term. Using p = mc'j(3 and the 
classical to quantum correspondence p —iKV we find: 

H ^ -ficr ■ \b - i{j+l)-^ — E xv\ , (25) 
I mc I 
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where we have parameterized the quantum evolution by 
the proper time, r. Thus, even at the classical field level, 
if one was to consider a spread in the momentum of the 
electron the Thomas term would alter the energy 

separation found in equation Furthermore, the electric 
field vacuum fluctuations perpendicular to the momentum 
would cause the particle to fluctuate about the classical path. 
Following the discussions relating to the perturbations of 
the classical path made in reference Jst] we assume, as they 
have, that the particle path remains classical and unperturbed. 
This seems to be a good approximation for particles with 
a gyromagnetic factor g ~ 2, such as electrons. Since we 
neglect perturbations to the path, the Thomas term in equation 
(I25I 1 should be dropped and the moving electron will remain 
within this approximation in a momentum eigenstate along 
the z-direction. 

Thus far we have only been concerned with the dynam- 
ics between ei and the electromagnetic fields, we now in- 
troduce the ancillary electron 62 into the Hamiltonian. We 
take the Hilbert space of the combined system to be 7Y — 
Ti-i (8) 7^2 ® ''Afield, where 7i 1 is the subspace of the accelerated 
spinor, 7^2 is the subspace of the spectator electron, and Tifieid 
is the Hilbert space for the electromagnetic field. In terms of 
the proper time the spin-field interaction Hamiltonian for the 
whole system is: 



H = -^la ®1- B. 



(26) 



Since the action on 7^2 is always the identity, for notational 
simplicity we suppress the factor of (Xil from all operators 
acting on the combined spin Hilbert space in what follows. 

We now consider the second quantisation of the magnetic 
field. We do this in the standard way by interpreting the mag- 
netic field as a space-time dependent field operator Expand- 
ing equation (l26T l. and defining the field operators, B± ~ 
Bx ± iBy, 1 19] we obtain the interaction in terms of field 
mode operators: 

V{x) = -n{(7xBx{x) + (TyBy{x) + (T^B^(x)), 

= -fiia^B+ix) + a+B^{x) + a,B,ix)), (27) 

where in the last step we have used a± = hjf^ x T 'i'^y)- Thus, 
the spin flips when a magnetic field mode 12011 is excited. 

Next we calculate the time evolution of the two electron 
system using a perturbative master equation ll2lll22|] whereby 
we interpret the magnetic field fluctuations about this system 
as an external environment ll23ll . As we are viewing the 
fluctuations of the external magnetic field as an unobserved 
environment of the spin-spin system, tracing over the field 
subspace is implied and thus no ad hoc dynamical modelling 
of the source of decoherence is necessary. 

We can write the total Hamiltonian of the system as 

H = Bsc + Bb + V, (28) 

where Hsc is the semiclassical evolution of the spin system, 
(fT9] l. Bb is the free Hamiltonian for the magnetic field. 



whose exact form will not be required, and V contains the 
quantum fluctuations arising from the second quantization of 
the spin-field interaction determined in equation ( |27] l. 

Let pt be the total density operator of the system plus field 
in the interaction picture. The equation of motion is: 



dT 



in 



(29) 



where r is the proper time of the accelerating electron. We 
assume that initially the electron spin system and field states 
are uncorrected so that: 



Pt{0) = p/(0) ® pB, 



(30) 



where pj describes the state of the two electron spins and pB 
describes the state of the magnetic field. Following [iT], one 
can expand equation (|29|l in a perturbative series and noting 
that T'<:b{V{t)pb) = 02411 one finds to second order in per- 
turbation theory: 



dpijr 

dT 



- / dnTVB[V{T),[V{T,),pi{T)®pB]]. 



For brevity we write the interaction in ( |27] | as 

V{T)^-pY.a,{T)Bl{x{r)), 



(31) 
(32) 



where i — {+, — , z} and the ai{T) are Heisenberg operators. 
By taking out the dependence of Bsc from pi using: 



P,(t) = e'"^^-/^ps{T)e 



'iHscT /h 



(33) 



we can rewrite the master equation (ISTT i for the spin subsys- 
tem in terms of ps, i.e., the reduced density matrix in the 
Schrodinger picture (see ll22ll '): 



dps{T) 

dT 



-[Hsc:Ps\ 



dri X 



TTB[a,Bj{x), [a,{Ti - t)b}{xi), ps ® pb]1 (34) 

where x = x{t) and xi = x{ti) are defined for short. This 
equation is still not entirely in the spin subsystem Schrodinger 
picture as we still have the aiji — r) Heisenberg operators. 
However since Bsc is given by equation ( fT9l ) we can use the 
Heisenberg equations of motion to find 



crj(r) 



(35) 



where aj — +1,-1,0, for j = +,— ,z, respectively, and 
defining the expectation value over the field state to be (O) = 
TrB(OpB) we obtain: 



dps 1 r„ p^ ^ 



+ {ps^ja, - (j,psaj){B]{xi)Bl{x))'^ . (36) 
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In Gaussian units the electromagnetic free field Wightman where j £ {x,y,z} ll26ll and 
function is 



\F^Ax)Fpa{x')\0) = {x - x')-'x 

TT 

-2 [(a; - x')p{x ~ x')pg^„ - {x ~ x')u{x ~ x')pg^^ 



-{x - x')f_,{x - x')„g^p + (a; - x!)v(x - x')a5Mp]} 



G{x,x') = —{x-x'Y^. 

TT 



(39) 



Along the constantly accelerated path described by equations 



(37-) (fT4li-(fT5]), this further reduces to: 



where F^^y — dpA^ — d^Ap is the anti-symmetric electromag- 
netic field strength tensor. Using = {(j>, A) and B = WxA 
one finds B = (F32, F13, F21). Since the motion is entirely 
along the z direction we find that the only non-zero two point 
correlation functions are 02511 



G(t-t') = G{x{t),x{t')) 



|sinh[^(T-r')]} . (40) 



{Bj{x)Bj{x')) — G{x{t ~ ie),x{T' j), (38) Thus we obtain the equation: 
I 



dps 1 r„ 1 
-^--f:[Hsc,Ps]-j^ 



2{a-a+ps - cr+pscr-) f dTie'^^''^-^^'''G{T ^ n - ie) + 2{psa+u- - a-ps(J+) f dTie'^^^'-^^/''G{T - n + ie) 
Jo Jo 

+ 2{a+a-ps-(T-ps(T+) f dTie-'^^''^-^^'^'G{T - n ~ ie) + 2{pscr-cr+ - a+ps(7-) f dTie''^'^^^-''^/^G{T - n 
Jo Jo 



+ (ps - <^zPs<^z) I {G{t - ti - if) + G{t - ti + ie)) dri 
Jo 



(41) 



Since the integrands in the above equation are sharply peaked where we have made use of the identity G{—z) — G{z) for 

functions about ri = r we are justified in making the Marko- any complex z, which follows from equation ( l40l i. This mas- 

vian approximation (see |22], pg.28). We make the change of ter equation is similar in form to those found in [22, 27J. By 

variables = t— ri and extend the integration over the interval separating Fi and r2 into their real and imaginary compo- 



s £ [0, 00]. To simplify the notation we define the following 
integrals: 



r± = 



G{s ± ie)ds, (42) 



{G{s-ie) + G{s + ie)}ds. (43) 



Using 



sinh*(s + ia) 



i coth(s + ia) (csch^(s + ia) — 2) , 



(44) 



we find that = fffB"- Then equation (Ell becomes: 



dps 1 r„ 1 r 
-^--^[Hsc,Ps]-j^{ 

2{a-a+ps - (T+ps(T-)T*^ + 2{ps<J+(T- - a^ps(T+)T 
+2{a+a^ps - cr-ps0'+)r_ + 2{pscr_a+ - a+ps(T-)T 
2ha^ 

3 TTC' 



+ o—iPs - crzPs(^z) }, 



nents equation ( |45T l simplifies into: 

dps 1 r„ 1 

-^ = -^[Hsc,Ps]~j^{ 

2Rer_(0-_(7+/9s + ps(7-(T+ - 2(T+ps(7-) 

+2Rer+((7+(T_ps + Pscr+cr- - 2ct_pso-+) 
2ha^ 

3 TTc' 



g [ps - cTzPs(^z) - «im(r+ + r_)[(T2, ps] } . 



(46) 



The last term is a correction to the unperturbed energy sepa- 
ration A. It is due to the spin-field coupling and effectively 

2 

renormalizes the Hamiltonian by a term H^c = h^^^ 



T_)az, reminiscent of the Lamb shift Il28ll in atomic physics. 
Further defining: 



7± ^ ^Rer, 



G{s-ie)ds, (47) 



(45) observe that the Master equation takes a manifestly Lind- 
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blad form: 



, (48) 



where {x, y} = xy + yx denotes an anticommutator, H' = 
Hsc — Hlc is the coherent part of the renormahzed spin 
Hamihonian, and Lj are the Lindblad operators, given by: 



L2 



2 

7+ 



(49) 
(50) 



(a pure dephasing channel) where 
,2 



By renormaHzing the subsystem density matrix. 



_ tH'r/h 



ps{T)e 



(51) 



(52) 



(53) 



we can write 



dp ^_ ^ ~ ~ 1 

+ ^ [2(T_|_pO-_ - Cr_CT_|_p — /5CT_tT+] 



(54) 



We note that this master equation is similar to the master 
equation for spontaneous emission of an atom discussed in 
129] (pg. 388) except that in our case we have a 4 x 4 density 
matrix and two more Lindblad operators. Nevertheless, 
we will now show that the same method that was used to 
solve the master equation in [29] can be applied here by 
choosing the right generalisation of the Bloch sphere in 4 x 4 
dimensions. 

Since p spans a sixteen dimensional vector space and the 
direct product of Pauli matrices including the identity, {cTi ^ 
(Jj\i,j £ 0, • • • , 3}, form sixteen linearly independent vec- 
tors we can expand any general density matrix for a two spin 
system as follows: 



3 3 

p^^^nja,(giaj, 
where we have chosen the Pauli matrices: 



(55) 



















I) 









1 

-1 

(56) 

and defined do = 1- A nice property about this choice of ba- 
sis is that the expansion coefficients are real, which follows 



from the hermiticity of the Pauli matrices and density opera- 
tor, furthermore the expansion coefficients can be computed 
directly using: 



1 



(57) 



As every density matrix has trace one, roo is equal to one quar- 
ter. From the inequality Tip^ < 1 the density operator can 
also be expressed as: 



p = ^ j l4x4 + ^ yijai 




for transitions down and up (in spin energy) respectively and where 



(58) 



(59) 



where equality holds if and only if the state is pure. Equation 
( |58] l generalises the Bloch sphere representation of a single 
qubit ll29ll to a two qubit system, in this way two qubit mixed 
states can be thought of as lying somewhere within a fifteen 
dimensional unit- sphere. 

In our case we will find it more useful to express our density 
in the form ( ISST l. As an example, a maximally entangled Bell 
state of the kind (| ft) + | H)) would be expressed 

PBcll j{cro CTo + 0-1 (g) (Ti - 0-2 0-2 + 0-3 (T3). (60) 

The entanglement between the spin degrees of freedom can 
be obtained by calculating the concurrence |01 



C{p) = max{Ai - A2 - A3 - A4, 0}, 



(61) 



where {Ai, A2, A3, A4} are the non-negative square roots of 
the eigenvalues of the matrix 



M = p{a2 (8) (y2)p*{a2 ® 02)- 



(62) 



Thus, the maximally entangled state (|60] l has a concurrence 
of one. 



Substituting p = j ^ij^i ® into equation ( |54] l gives: 



dr. 



11 7- 
— — fjj (K) CTj = — [2cr_crjCT+ - (T+o-_cri - criCr+cr_J 



, 7+~ 



(63) 



which after a little algebra gives sixteen first order linear dif- 
ferential equations 

roj{r) = 0, 

hjir) = -i(7_ +7+ +472)fij(r), 

hjir) = -i(7_ +7+ +472)r2i(r), 

hj{r) = (7- - 7+)^0j(t) - (7- +7+)f3j(T), (64) 
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where dots imply differentiation with respect to t. The solu- 
tions to these equations are found to be: 

rojir) = ?0j(0), 



^r~o,(0)(l-e-(^-+-'+)^). 



(65) 



Thus, the relaxation and dephasing times respectively are: 



t; 



7- +7+ = 
7- +7+ 



+ 27., 



(66) 
(67) 



where T2 is the timescale for which the off-diagonal ele- 
ments of the density matrix ('coherences') decay and Ti is 
the timescale for spin flipping. The relations (|66]|-(|67]| have 
a form identical to those typically found in magnetic reso- 
nance problems (see |27], chapter 4. 3). We will compare these 
timescales in more detail in the next section once we have cal- 
culated the spin flip rates. 



IV. RELAXATION AND CONCURRENCE 

We have succeeded in finding the time evolution for the 
system of spins within our specific setup. In order to analyse 
the effect that the acceleration has on the entanglement shared 
by these spins it is necessary to calculate the spin flip rates 
defined in equation ( |47] i. 

Making a linear change of variables, s' — as/c, the integral 
in ( |47] ) becomes 



7± 



2/i2 ha'^ c 



°° ,exp (=FzAs'c/a/i) 
as 



sinh*^ i(s' — ie) 



(68) 



These integrals can be evaluated using contour integration Jsj] : 



7± = To 



Bell and Leinaas also noted that the ratio of the transition 
rates, 



7+ _ -2TrcA/ah 

7- 



(70) 



define an equilibrium ratio of populations of the upper and 
lower states. Thus, the equilibrium distribution over the levels 
has a thermal character in accordance with the Unruh temper- 
ature formula ([T]i- To proceed further we recall that the Bose 
occupation number is 



1 

o2-!TcA/ah _ ^ ■ 



(71) 




FIG. 2: The relaxation times Ti (dashed line) and T2 (solid line) 
in units of the spontaneous emission timescale at zero acceleration, 
7(^^, as a function of the acceleration, a, in units of cA/fi. 



We also notice from equation f69i that there remains an emis- 
sion rate even when the acceleration goes to zero: 



70 = lim7_ = ij-r^. 

a^o 3 n^c-' 



(72) 



This spontaneous emission is due to the quantum interactions 
of the magnetic moment with the magnetic field and defines a 
useful parameter, 70, which sets the scale of spin flip transi- 
tion rates. This magnetic spontaneous emission rate is much 
weaker than the usual atomic spontaneous emission rate, and 
is significantly harder to detect, see I,3Q1 . The equations in 
(l47T i can now be written: 



7+ 



70 



n, 



7- = 7o(l + (^)')(n+l). 
Furthermore, it is easily verified that 



(73) 
(74) 



1z — 4^ 70 ( . 



(75) 



Thus, the relaxation and dephasing times respectively are: 



1 - 7o(l + (ffr)coth^, (76) 
Ti' = ^{(l + (^)^)coth + 1(^)3}. (77) 

The relaxation and dephasing times are shown in figure|2] We 
observe that Ti < T2 < 2Ti for all values of the acceleration. 
Thus, typically there will be more than one spin flip before the 
coherences vanish. 

In the long time limit the state ( fSSl l becomes: 

Pa(oo) = ^fojl«)o-j+tanh(^)fojCr3«)CTj, (78) 

j 

(l + tanh(^)a3) 



(79) 
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which is a product of two density matrices acting separately 
on each of the spinor subspaces. Therefore the state of the 
spins will eventually separate regardless of how they were 
initialised. To determine precisely how the system separates 
we calculate the concurrence (|6T1|-(|62]| as a function of time. 
Before doing so, it is worth discussing the effect of Lorentz 
transformations on entanglement. 

Recall that Gingrich and Adami 1^ have shown that the 
concurrence of two particles, each in a state with some spread 
in momentum, is not Lorentz invariant. This is because 
momentum dependent Wigner rotations act on the spinors 
under Lorentz transformations (see also 1 18] for a discussion 
on the entropy of single particles under Lorentz transforma- 
tions). One then wonders if in the present case calculating 
the concurrence in the rest frame of the inertial electron is 
meaningful. While this is clearly a concern for particles 
with a spread of momentum states, when the particles are in 
momentum eigenstates the Wigner rotation acts like a local 
unitary operation on the spinors and therefore does not change 
the concurrence ITit]. In our case, the moving electron 
follows the classical path defined by equations (fT4b-(fT5l) and 
is thus by construction in a momentum eigenstate Pz{t)- 
Furthermore, the stationary spinor can be chosen to be as 
narrow as is required i.e., ~ 0. Under these conditions, 
since there is no significant momentum spread in either 
of the particles wavefunctions, the concurrence function 
behaves to good approximation invariantly under Lorentz 
transformations. We now proceed to calculate the value of 
the concurrence in the stationary electron rest frame. 

We initialise the system into the maximally entangled Bell 
state ( |60l l which then evolves according to equation ( |65] l. since 
concurrence is invariant under local unitary transformations 
|j7|] we are justified in calculating the concurrence of p instead 
of PS- The time-dependent system density matrix is: 



~P{r) 



{(To (Xi (To + e '"^^CTi 



I (Tl — e ^^'^ (72 



(72 



-e-^'^73 



(T3 



\ ah J 



e ^1^) erg® (To} 



(80) 



To"\ Oneob- 



where we have defined Fi = Tj~^ and T2 — -L2 
serves from equation dSOl ) that p(t) = p*{t), and since the 
eigenvalues of M in equation ( |62] | are real and positive |7] we 
can find the concurrence of /5(t) by diagonalising the matrix: 



pa2 (XI (72 , 



(81) 



where we take the {Ai} in equation dMT l to be the absolute 
values of the eigenvalues of the matrix dSTT i. We find that the 
concurrence is given by: 



C(a, 



{ 



1 - e 



-r^l^ 



sech(^),0} 



(82) 

The concurrence is plotted in figure [3] as a function of the ac- 
celeration and the proper time. We observe that the greater 
the acceleration the quicker that the initial entanglement dis- 
appears. We will use the time taken for the system to reach 
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FIG. 3: (Color online) We have plotted the concurrence of two ini- 
tially maximally entangled spins as a function of the acceleration, a, 
(in units of ^) and the proper time (in units of Jq^). Our main re- 
sult, equation ([85}, is also shown overlaid in a dark (red) line along 
C = 0. 



zero concurrence to quantify the dependence of disentangle- 
ment on acceleration. For a given acceleration the time taken. 
To, for the system to completely disentangle is given by the 
equation: 

g-ror2 = i (1 - e-^"^^) sech(i^), (83) 

which follows from equation (|82l ) jsill . Defining the dimen- 
sionless parameter a = ^ and taking the limit a ^ 1 (i.e., 
A < ^) we find, 

8 p'^a^ 



Fi = Fa = 



Equation 



can now be solved for tq: 

37rln3 ftc^ 



0{a). 



To 



(84) 



(85) 



Thus in the large acceleration small magnetic field limit the 
proper time taken to disentangle the two spinors is propor- 
tional to the inverse of the acceleration cubed [32]. We have 
used this result to plot the zero concurrence dark (red) line in 
figure [3] 

In the frame of the stationary spinor the time taken to disen- 
tangle is exponentially longer: 

'SvrlnS hc^ 



tn 



8 



p^o? 



If we put in the numbers for an electron we find: 



la 



3.8 X lO^^^s 



(86) 



(87) 



Therefore an acceleration with a magnitude of about thirty 
is required to observe the disentanglement on a reasonable 
timescale. This is consistent with the thermal equilibrium 
timescale found in jsj]. 
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V. CONCLUSION 

We have considered the system of two spin-entangled elec- 
trons when one electron is accelerated and placed under a con- 
stant magnetic field whilst the other is at rest and isolated. 
Our method consisted of explicitly calculating the open quan- 
tum system where the quantised magnetic field fluctuations 
were considered to be an unobserved environment. This gen- 
eralises the linearly accelerated single electron case consid- 
ered by Bell and Leinaas over 25 years ago. For the first time 
we have found analytic expressions for the T2 dTTl l. the con- 
currence ( l82b and the entanglement lifetime ( |86] | of this sys- 
tem. 

It is worth emphasising that while the timescales found here 
for the linear system prohibit any experimentation, there is a 
drastic improvement in the timescale when the accelerating 
electron is put into a circular orbit |5]. Using the methods we 
have employed here, one could perform the calculation in the 
case when the motion is circular or when the acceleration in 
the IRF is simple harmonic. The latter situation is of interest 
as it would be a good theoretical basis in which to study the 
entanglement between spins attached to high frequency can- 
tilevers. 



Acknowledgments 



Using R[dn)A{d(3') = A{dj3')R{dn) + 0{dl3) S one 
finds: 

R{dn) = A{-d0')A{0 + d0)A-^{0), (A3) 



= I-dn-S, 
where S is the generator of rotations and 



dn 



7 



7 + 1 



/3 X dl3. 



(A4) 



(A5) 



One can verify that —d/S' is the three component vector of 
A{(3 + d(3)f3, where (3 = (7, 7/3), in what follows we adopt 
the notation that (3\ represents the three vector component of 
A/3. Defining L = A^^ and A = A(0 + d0) we have. 



R{dn)^L-\0tOKL{0), 



(A6) 



which shows that the infinitesimal Thomas rotation is a special 
case of the Wigner rotation f^Jl^ for the choice of A given 
above. Therefore a state \p\) acted on by U{K) (where U 
is a two dimensional linear operator satisfying [/(A1A2) = 
C/(Ai)[/(A2)) can be written 
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= Y,Dxx'{K,p)\pk\'). (A9) 
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APPENDIX A: CONNECTION BETWEEN THOMAS 
PRECESSION AND THE INFINITESIMAL WIGNER 
ROTATION 



As discussed in the introduction a single boost A{l3 + d(3) 
differs from the combination of boosts A{(3) followed by 
A{d(3') by a rotation: 

R{dVl) = A{0 + d0)[A{d0')A{0)]-^ , (Al) 

where rf/3' is measured in the IRF and is related to d/3 = v /c 
(see LMl, chapter 11.8, LMJ) by: 



d0' = jd0- 



r 



7 + 1 



d(3 ■ (313. 



(A2) 



where we have made use of L{f3)'^k'^ — p^, where k is the rest 
frame four momentum (toc^, 0). On the second line we have 
inserted the identity U{L{0^))U{L~'^{0k)) 1 and identi- 
fied the Wigner Rotation. The momentum dependent Wigner 
D-function is known to disentangle spin states iQl when ei- 
ther state has some spread in momentum. Putting the rotation 
into this form allows us to see that a continuous sequence of 
Wigner rotations caused by a constant acceleration in general 
leads to a mixing of the spin and momentum entanglement. 
We now complete the analysis by showing that this effect is 
just the Thomas precession, which can be accounted for by 
including an additional term in the Hamiltonian. 
When the acceleration in the IRF is constant the moving frame 
rotates with a constant angular velocity. 



Now we can write 



dfl _ 72 d/3 
dt ^ j+1 dt 

dn 



X /3. 



U(R(dn)) = 1 - — . m, 
dt 



(AlO) 



(All) 



where in our case we choose § to be the generators of the 
two dimensional representations of the unitary subgroup of 
5L2(C ), i.e., that satisfy the commutation relations: 



^ijk^k, 



(A12) 
(A13) 
(A14) 
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where K are the hermitian subset of S'L2(C ) corresponding 
to boosts. Since the Pauli matrices satisfy [cTj, a-j] = 2zeyfe(Tfe 
we can take, 



Then 



-, yie{x,y,z}. 



U{R{dn)) = 1+ujT-S, 

i ( hujT • 



= 1 - - 



h\ 2 



(A15) 

(A16) 
dt. (A17) 



This equation is the infinitesimal form of the unitary time 
translation operator and thus the generator in brackets is the 
Thomas Hamiltonian, 



hujT ■ <y 



(A18) 



In the IRF the spin-field equation of motion is given by. 



^= flxB' (A19) 
or 



where B' is the field in the IRF. In the laboratory time. 



^ = 7- V X B\ (A20) 



Thus, the total Hamiltonian for the relativistic two component 
spinor can then be written. 



H = -7"V- ^ + — 7; • 



(A21) 



[1] S. W. Hawking. Commun. Math. Phys., 43:199-220, 1975. 
[2] P. C. W. Davies. /. Phys., A8:609-616, 1975. 
[3] Stephen A. Fulling. Phys. Rev., D7:2850~2862, 1973. 
[4] W. G. Unruh. Phys. Rev., D14:870, 1976. 
[5] J. S. Bell and J. M. Leinaas. Nucl. Phys., B212:131, 1983. 
[6] Emil T. Akhmedov and Douglas Singleton. 2006. hep- 
ph/06 10391 

[7] Wootters,W.K.. Phys. Rev Lett., 80:2245-2248, 1998. 

[8] Robert M. Gingrich and Christoph Adami. Phys. Rev. Lett., 

89(27):270402, Dec 2002. 
[9] Paul M. Alsing, David McMahon, and G. J. Milbum. J. Opt. 

B: Quantum Semiclass. Opt. 6 (2004) S834-S843, arXiv:quant- 

ph/0311096 vl 

[10] Alsing,PM. Milbum,G.J. Phys. Rev Lett, 91:180404, 2003. 
[11] P. M. Alsing, Ivette Fuentes-SchuUer, Robert B. Mann, and 

T. E. Tessier. Phys. Lett., A74:032326, 2006. 
[12] Jian-Ming Cai, Zheng- Wei Zhou, Ye-Fei Yuan, and Guang-Can 

Guo. Phys. Rev A, 76:042101, 2007. 
[13] Fuentes-Schuller,!. Mann,R.B. Phys. Rev A, 95:120404, 2005. 
[14] Alsing, P.M. and Milbum, G. J. Quantum Inf. Comput. 2, 487 

(2002). 

[15] J. D. Jackson. Rev Mod. Phys., 48:417-433, 1976. 
[16] J. D. Jackson. Classical Electrodynamics Third Edition,Wiley, 
New York, 1998. 

[17] In [5] this effective Hamiltonian is written in terms of the IRF 
electric field. One can show using the Lorentz force equation 
(under the assumption that non-electromagnetic causes of the 
acceleration are negligible), that these two forms are equivalent. 

[18] Peres, Asher and Scudo, Petra F. and Terno, Daniel R. Phys. 
^ev. Lerf., 88, 23040, 2002. 

[19] This is just a change of basis of the standard polarization direc- 
tions and should not be confused with the creation or annihi- 
lation operators which are obtained by separating the field into 
positive and negative frequencies. 

[20] These magnetic modes are related to photon modes through 
B = \/ y. A. 

[21] Walls D.F. and Milbum G.J. Quantum Optics. Springer Verlag, 
Berlin, 1994. 



[22] Juan Pablo Paz and Wojciech Hubert Zurek. quant-plt/OOlOOl 1 , 
2000. 

[23] Since the spin-field coupling strength is small a perturbative ex- 
pansion is acceptable. 

[24] 1rB{V{T)pB) = -mE '^»Trfl(Bjpi3) and since B, = B+ + 
B~ where Bf are the positive and negative frequencies, then 
as long as the field is not in a squeezed state (we assume ps is 
not)Trfl(B±pi3) =0. 

[25] Where the choice of — ie has been made to recover the appro- 
priate Greens function for this two-point correlation function 
(see issll for a discussion on the analytic properties of Greens 
functions). 

[26] Note that {B±{x)B^{x')) = {B^{x)B^{x')) + 

(By{x)By{x')) = 2G{x,x'). 
[27] A.J. Fisher. Graduate course on open quantum systems. 2004. 

|www . cmmp . ucl ■ ac . uk/""aj f /course-notes ■ pdf | 
[28] J. J. Sakarai. Advanced Quantum Mechanics. Addision Wesley, 

1967. 

[29] M. A. Nielsen and I. L. Chuang. Quantum Computation and 
Quantum Information. Cambridge, 2000. 

[30] Gover,A. Phys. Rev Lett., 96(12): 124801, 2006. 

[31] In practice one could observe a loss of concurrence before this 
time, however, defining the timescale in this way is preferred 
since it is more meaningful to speak of the time taken for the 
spins to completely disentangle from each other as opposed to 
the time taken for them to evolve to some arbitrary value of 
concurrence. 

[32] We note for completeness that the time taken for the con- 
currence to decay to the value C, is given by rc = 

8 \2C+1 ) JiT^- 

[33] N. D. Birrell, P. C. W. Davies Quantum Fields in Curved Space 
Cambridge 

[34] However notice the slight difference in notation, df3' A/3. 

[35] See equation (11.117) of [16] although note the typographi- 
cal error in the argument of the last boost which should read 
Ai,oost(A/3) instead of Aboo3t{S(3). This error has been con- 
firmed via private correspondence with J. D. Jackson. 



